Environmental Noise and Nonlinear 
Relaxation in Biological Systems 

B. Spagnolo^^ D. Valenti^} S. Spezia}}L. Curcio^, 
N. Pizzolato^} A. A. Dubkov^^A. Fiasconaro^'^, D. Persano Adorno^, 
P. Lo Bue^, E. Pert" and S. Colazza^ 
^ Dipartimento di Fisica, 
Group of Interdisciplinary Physics 
Universita di Palermo and CNISM-INFM, Unita di Palermo 
Viale delle Scienze, ed. 18, 1-90128 Palermo, Italy 
^ Dipartimento di Ingegneria Elettrica, Elettronica e delle Telecomunicazioni, 
Viale delle Scienze, ed. 9, 1-90128 Palermo, Italy 
^ Radiophysics Department, Nizhniy Novgorod State University, 
23 Gagarin ave., 603950 Nizhniy Novgorod, Russia 
Departamento de Ffsica de la Materia Condensada, Universidad de Zaragoza, 

E-50009 Zaragoza, Spain 
^ Dipartimento di Scienze Entomologiche, Fitopatologiche, 
Microbiologiche, Agrarie e Zootecniche, Universita di Palermo, 
Viale delle Scienze, ed. 5, 1-90128 Palermo, Italy 



Abstract 

We analyse the effects of environmental noise in three different biological sys- 
tems: (i) mating behaviour of individuals of Nezara viridula (L.) (Heteroptera Pentato- 
midae); (ii) polymer translocation in crowded solution; (iii) an ecosystem described 
by a Verhulst model with a multiplicative Levy noise. Specifically, we report on 
experiments on the behavioural response of N. viridula individuals to sub-threshold 
deterministic signals in the presence of noise. We analyse the insect response by 
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directionality tests performed on a group of male individuals at different noise in- 
tensities. The percentage of insects which react to the sub-threshold signal shows a 
non-monotonic behavior, characterized by the presence of a maximum, for increasing 
values of the noise intensity. This is the signature of the non-dynamical stochastic res- 
onance phenomenon. By using a "hard" threshold model we find that the maximum of 
the signal-to-noise ratio occurs in the same range of noise intensity values for which 
the behavioral activation shows a maximum. In the second system, the noise driven 
translocation of short polymers in crowded solutions is analyzed. An improved version 
of the Rouse model for a flexible polymer has been adopted to mimic the molecular 
dynamics, by taking into account both the interactions between adjacent monomers 
and introducing a Lennard-Jones potential between non-adjacent beads. A bending 
recoil torque has also been included in our model. The polymer dynamics is simu- 
lated in a two-dimensional domain by numerically solving the Langevin equations of 
motion. Thermal fluctuations are taken into account by introducing a Gaussian un- 
corrected noise. The mean first translocation time of the polymer centre of inertia 
shows a minimum as a function of the frequency of the oscillating forcing field. In the 
third ecosystem, the transient dynamics of the Verhulst model perturbed by arbitrary 
non-Gaussian white noise is investigated. Based on the infinitely divisible distribution 
of the Levy process we study the nonlinear relaxation of the population density for 
three cases of white non-Gaussian noise: (i) shot noise, (ii) noise with a probabiUty 
density of increments expressed in terms of Gamma function, and (iii) Cauchy stable 
noise. We obtain exact results for the probability distribution of the population density 
in all cases, and for Cauchy stable noise the exact expression of the nonlinear relax- 
ation time is derived. Moreover starting from an initial delta function distribution, we 
find a transition induced by the multiplicative Levy noise, from a trimodal probabil- 
ity distribution to a bimodal probability distribution in asymptotics. Finally we find a 
nonmonotonic behavior of the nonlinear relaxation time as a function of the Cauchy 
stable noise intensity. 

PACS: 87.18.Tt, 87.50.yg, 05.40.-a,64.70.km 

1. Introduction 

During last decades noise-induced effects have been experimentally observed and theoreti- 
cally studied in different physical and biological contexts IIllHSIlllSllHlTlllim, such as 
neuronal cells, excitable systems and threshold physical systems |[T0l[TTl[T2l[T3l[T4l[T5llT6l 
[I7l[l8l[lll20ll2ll- 

In particular, stochastic resonance, resonant activation and noise enhanced stability phe- 
nomena in neuronal activation have been recently discussed 11211 l22l|23l . 

Nature consists of open systems characterized by intrinsically non-linear interactions 
and subject to environmental noise ll24l . The presence of random fluctuations, that are an 
uneliminable component of natural ecosystems, makes difficult detection and transmission 
of signals and can modify the information transported. 

However, in the presence of some specific non-linearity of the system and for suitable 
intensity of noise, counterintuitive phenomena, such as stochastic resonance (SR), can be 
observed. This indicates that noise can play a constructive role, improving the conditions 
for signal detection. 

SR phenomenon initially was observed in the temperature cycles of the Earth ll25l . can 
be found in many physical and biological non-linear systems |[26ll27ll28l . SR can be mod- 
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elled by a bistable potential subject to periodical driving force in the presence of external 
additive noise. The signature of SR is a non-monotonic behaviour, characterized by a max- 
imum, of the signal-to-noise (SNR) ratio as a function of the noise intensity. This indicates 
that the noise can enhance the amplitude of deterministic signals, improving the response of 
the system thi'ough a resonance-like phenomenon lfTTl[T2l[T3l[T5l[M[r7l[T8l[T9llMl2Tl[29l 
[30l[ni[3l[32[Ml[35l[36l|2i|27l[37l- However, SR does not occur only in bistable systems, 
but also in monostable, excitable, and non-dynamical systems. In these situations we name 
this effect non-dynamical (or threshold) stochastic resonance, because the phenomenon is 
connected with the crossing of a threshold and can occur also in the absence of an exter- 
nal potential HB [El |29l. Sensory neurons, that are threshold systems characterized by 
intrinsic noise, are an ideal workbench to observe non-dynamical SR because ll30l [3111321 . 
Historical experiments revealed the presence of non-dynamical SR in the neural response 
of mechanoreceptor cells of crayfish |[33l . and the improvement of sensorial activity of 
paddlefish in the detection of electric signals produced by preys |[34l l35l [36l . Such sen- 
sory neurons are ideally suited to exhibit SR as they are intrinsically noisy and operate as 
threshold systems ll30l[3n[32l . 

In this contribution, we study the effects of external noise in three different biological 
systems. We stait analyzing the mating behaviour of individuals of A^. viridula (L.) (Het- 
eroptera Pentatomidae). In particular, we investigate the role of noise in the response of 
male insects to mechanical vibrations emitted by female individuals and transmitted in the 
substrate |[38l l39l l40l . A'^. viridula, the southern green stink bug, is a pentatomid insect 
highly polyphagous and quite harmful for agricolture PTl l42l . N. viridula has up to five 
generations per year Il43ll44ll45ll46l . 

The mating behavior of A^. viridula can be divided into long-range location and short- 
range courtship. The first one includes those components of the behavior that lead to the 
arrival of females in the vicinity of males. The long range attraction mediated by male 
attractant pheromone enables both sexes to reach the same plant. 

Here, we analyze the mating behaviour of insects during the short-range courtship, 
when bugs of both sexes are very close and the acoustic stimuli (improperly called songs) 
can be an important element in the sexual communication ll38l . 

The sound is produced by the tymbal, an organ sited in the back and present in adult 
individuals |[39l . The vibrations, produced by a bug at the frequency of about 100 Hz, prop- 
agate through the legs into the plant stem and can be detected by the vibro-receptors placed 
in the legs of another insect B31I471 . Many experimental studies have been performed on 
this acoustic communication, analysing the different signals characteristic of populations of 
N. viridula from Slovenia, Florida, Japan and Australia BSl . 

The fundamental role of the vibratory signals suggests that a better knowledge of the 
mechanism of acoustic communication during the short-range courtship can help to point 
out more efficient strategy to control A'^. viridula populations, devising "biologic" traps 
whose working principle is the emission of acoustic signals. In natural conditions, A'^. 
viridula populations interact strongly with environment, and therefore the presence of sur- 
rounding noise becomes an essential component of the acoustic communication. 

In the second part of this contribution, we consider transport phenomena of polymers in 
crowded solutions. In fact, the translocation of DNA and RNA across nuclear membranes as 
well as the crossing of potential barriers by many proteins represents a fundamental process 
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in cellular biology. The study of the transport of macromolecules across nanometer size 
channels is important for both medical research in anticancer targeted therapy B9l l50l and 
technological applications |[5Tl[52l . 

First experiments on the passage of DNA molecules across an a-hemolysin (a-HL) 
protein channel revealed a linear- relationship of the most probable crossing time Tp with 
the molecule length ||53l . Moreover, Tp scales as the inverse square of the temperature and 
the dynamics of biopolymer translocation across an a-HL channel is found to be governed 
by pore-molecule interactions Il54ll55l [56l. More recent experimental studies have shown 
that the application of an AC voltage to drive the translocation process of DNA molecules 
through a nanopore plays a significant role in the DNA-nanopore interaction, and provides 
new insights into the DNA conformations |[57ll58ll59ll60ll6TI . 

The complex scenario of the translocation dynamics coming from experiments has been 
enriched by several theoretical and simulative studies ||62l|62|64ll65l|6ll|671|68l|69l. The 
mean first passage time of a Brownian particle to cross a potential barrier in the presence of 
thermal fluctuations and a periodic forcing field has been theoretically and experimentally 
investigated as a function of the driving frequency in Refs. 1170117111721 1731 1741 1751 . The 
translocation time of chain polymers has been theoretically studied in the presence of a di- 
chotomically fluctuating chemical potential only as a function of its amplitude in Ref. [76 1. 

In particular we investigate the role of an external oscillating forcing field on the trans- 
port dynamics of short polymers surmounting a barrier, in the presence of a metastable 
state. We find a minimum of the mean first translocation time (MFTT) of the molecule 
center of mass as a function of the frequency of the forcing field. This nonlinear behaviour 
represents the resonant activation (RA) phenomenon in polymer translocation. We find that 
a suitable tuned oscillating field can speed up or slow down the mean time of the translo- 
cation process of a molecule crossing a bamer, using the frequency as a control parameter. 
This effect can be of fundamental importance for all those experiments on cell metabolism, 
DNA-RNA sequencing and drug delivery mechanism in anti-cancer therapy. 

In the third part of this chapter we investigate the transient dynamics of the Verhulst 
model perturbed by arbitrary non-Gaussian white noise. The nonlinear stochastic systems 
with noise excitation have attracted extensive attention and the concept of noise-induced 
transitions has got a wide variety of applications in physics, chemistry, and biology [82]. 
Noise-induced transitions are conventionally defined in terms of changes in the number of 
extrema in the probability distribution of a system variable and may depend both quanti- 
tatively and qualitatively on the character of the noise, i.e. on the properties of stochastic 
process which describes the noise excitation. The Verhulst model, which is a cornerstone 
of empirical and theoretical ecology, is one of the classic examples of self-organization in 
many natural and artificial systems ll83l . This model, also known as the logistic model, is 
relevant to a wide range of situations including population dynamics Il82ll84ll85l[86l . self- 
replication of macromolecules ll87l . spread of viral epidemics [i88l . cancer cell population 
ll89l . biological and biochemical systems liOOllOll . population of photons in a single mode 
laser |[92l[93l . autocatalytic chemical reactions |[94l[95ll96ll97ll98l . freezing of supercooled 
liquids ||99l, social sciences IIIOOIIIOIII . etc. 

By considering the season fluctuations and the random availability of resources we an- 
alyze the stochastic Verhulst equation in the presence of a non-Gaussian stochastic process. 
By investigating the transient dynamics of this model we obtain exact results for the mean 
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value of the population density and its nonstationary probability distribution for different 
types of white non-Gaussian. Noise-induced transitions for the probability distribution of 
the population density and a nonmonotonic behavior of the nonlinear relaxation time as a 
function of the Cauchy noise intensity are found. 

The chapter is organized as follows. In section 12X1 we report on experimental setup and 
methods used in the investigation of behavioural response in A^. viridula. In section l2!2l we 
present our experimental results of directionality tests on the behaviour of male individuals 
of N. viridula. In section [231 we discuss the experimental findings and compare them with 
theoretical results obtained by a hard threshold model. 

In Sect. |3] we present our polymer chain model and give the details of the molecular 
dynamics simulation process. Results are reported in Sect. 13.1.1 In the next section IV 
we present our Verhulst stochastic model with Levy noise excitation together with all the 
theoretical results obtained. Finally conclusions are drawn in Sect. 171 

2. Behavioural Response in A^. viridula 
2.1. Materials and methods 

In our experiments we used individuals of N. viridula collected in the countryside around 
Palermo, and reared in laboratory conditions 1791 . Male insects have been used for experi- 
mental trials after they reached sexual maturity (not less than ten days after the final moult), 
and a thiee-day period of isolation from the opposite sex B0ll48l . 

The sexual calling song emitted from a female individual has been recorded by the 
membrane of a conic low-middle frequency loudspeaker (MONACOR SPH 165 C CAR- 
BON with a diameter of 16.5 cm). Afterwards the sound, stored on a pc, has been analysed 
and processed using a commercial software. The speaker has been used as an "inverse" 
microphone, namely an acoustic-electric transducer: the sounds have been recorded from a 
low-frequency non-resonating membrane of a speaker, conveniently chosen to get a good 
frequency response starting at 20 Hz. The sound acquisitions have been made inside an 
anechoic chamber (sound insulated) at 22 — 26°C, 70 — 80% of relative moisture and in 
presence of artificial light. The choice of this recording set-up has been decided after a 
comparative analysis with a recording system based on the use of a stethoscope. In partic- 
ular, the speaker membrane shows greater sensitivity at medium-low frequencies, that are 
crucial to our experiment. 

The sound has been sampled from the analogical signal source (44100 samples per 
second at 16-bit) and then filtered by an IS*'^ order Tchebychev filter (type I) with band- 
pass from 60 to 400 Hz. This filtering has been done to cut: (i) the low frequencies due to the 
electric network (50Hz) and those from the conic loudspeaker, and (ii) the high frequencies 
due to the electronic apparatus. Spectral and temporal properties of the measured non- 
pulsed female calling songs (NPFCS) have been compared with those of North America, 
observing that A^. viridula individuals collected in Sicily have the same dialect as adults of 
A^. viridula collected in USA with a slightly different frequency range B0ll48l[80l . 

In Fig. [It, the oscillogram of NPFCS is shown. The signal is characterized by a short 
pre-pulse followed by a longer one, according to previous experimental findings BSl . In 
Fig. Wp, the power spectrum density (PSD) of NPFCS is shown. In this spectrum the dom- 
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Figure 1. (a) Oscillogram (b) Power spectrum density and (c) Sonagram of the non pulsed 
type of Nezara viridula female calling song; (d) The block diagram of the experimental 
setup. 



inant frequencies range from 70 to 170 Hz and the subdominant peaks do not exceed 400 
Hz. The maximum peak occurs at 102.5 Hz. In Fig. [B; we report the relative sonagram, 
achieved by the Short Time Fourier Transform (STFT) method. The STFT maps a signal 
providing information both about frequencies and occun^ence times. It shows that during 
the first two seconds (short pre-pulse) the dominant frequency interval is narrower than the 
range observed in the subsequent time space. In particular in the first time interval the high- 
est frequency does not exceed 130 Hz, whereas in the final one it reaches almost 170 Hz. 

Here we study the effects of noise on the behavior of A^. viridula during the mating 
period. Therefore, in order to perform directionality tests, we have designed and constructed 
a Y-shaped dummy plant, and placed it inside an anechoic chamber. The Y-shaped plant 
consists of a vertical wood stem, 10 cm long, and 0.8 - 0.9 cm thick at the top of which there 
are two identical wooden branches, 25 cm long, and 0.4 cm thick, as shown in Fig. [T]l. The 
angle between two branches is 30° — 50°. 

In our experiment a signal is sent along one branch of the Y-shaped substrate and the 
behavior of single male individuals, initially placed at the center of the vertical stem, is 
observed ll38l (see Fig.[TJl). The source of vibratory signals (i.e. the cone used as an electro- 
acoustic transducer) is in contact with the right apex of the Y-shaped dummy plant. Vertical 
stem and lateral branches aie not in direct contact, albeit in close (0.5 cm) proximity. 

We consider a trial valid if the insect, before choosing one direction in the Y-shaped 
structure, has checked the two possible directions of signal origin, touching the lower ex- 
tremity of both branches. By following these criteria, we made several observations for 
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different intensities of female calling songs, recording the choices (left or right) of each 
male individual used in our trials, and obtaining a set of statistical data that allows us to 
determine the intensity threshold value at which the bugs start to "hear" the calling song. 

2.2. Experimental results 

The presence of an "oriented" behavior, that is the tendency of the insects to choose the 
branch with the signal source, is revealed by performing directionality tests on a group of 
male individuals. When we observe a percentage of insects higher than 65% going towards 
the acoustic source, source-direction movement (SDM), we consider that the signal has 
been revealed by the insects. In Fig. we plot the relative frequency of SDMs, that is 
the number of SDMs divided by the total trials, at different signal intensities. The exact 
number of trials, performed for each intensity, is reported beside the corresponding point in 
the graph. For small values (lower than 0.0010 V) of the signal power approximately 50% 
of the insects choose one direction and the remaining 50% the other. Conversely, for values 
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RMS Amplitude [V] ' ' 

Figure 2. Plots of the Source-Direction Movement (SDM) Ratio as a function of: (a) the 
female calling song Root Mean Square (RMS) amplitude (purely deterministic signal); (b) 
the noise intensity D. In each experimental value is reported the error bar and beside the 
corresponding number of the performed trials. 



greater than 0.0020 V, the insects show a preferential behavior, choosing the direction from 
which the signal originates in the 80% of the trials. Consequently, we have chosen the value 
0.0015 V of the signal power as the threshold level for signal detection. 

Then, by using a sub-threshold signal plus a Gaussian white noise signal we have in- 
vestigated the response of the test insect for different levels of noise intensity D. In Fig. [2}) 
we report the percentage of SDMs as a function of D, finding the optimal noise intensity 
that maximizes the recognition between individuals of opposite sex. The graph shows a 
maximum for D ^ 1.30 • 10"^ V^. For values of D both lower and higher than 1.30 • 10""^ 
V"^, the response of insects is not significant. In particular, for D > 1.30 • 10~^ V"^ the 
percentage of individuals going towards the acoustic source decreases below 0.5 reaching 
0.2 for D ^ 3.75 • 10"^ F^. The other values of the SDM ratio close to 50%, indicate that 
individuals of A^. viridula randomly choose the direction of their motion, that is no oriented 
behavior occurs. The non-monotonic behavior of SDM, with a maximum at D w 1.3010^^ 
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V'^, indicates that in the presence of a sub-threshold deterministic signal, the environmental 
noise can play a constructive role, amplifying the weak input signal and contributing to im- 
prove the communication among individuals of N. viridula. The occurrence of a minimum 
in the SDM behavior at L) a; 3.75 • 10~^ V"^, will be subject of further investigations. A 
possible conjectural explanation is the following: when the noise intensity is so large that 
the signal received from the vibro-receptors is significantly modified, the male insects are 
not able to recognize the female calling song, and they exchange it for the song of some 
rivals. 

A further increase of the noise intensity causes the spectrum of the received signal to 
become indistinguishable from a pure environmental noise and therefore the insect is unable 
to recognize any signal of A^. viridula individuals. This implies that no significant response 
is observed in terms of percentage of source-direction movements (SDMs ~ 50%). 

2.3. Threshold Stochastic Resonance 

The presence of a maximum in the behavior of SDM percentage as a function of D can 
be explained either by the threshold phenomenon, or non-dynamical, stochastic resonance 
(TSR). 

Stochastic resonance (SR), initially observed in the temperature cycles of the Earth 
ll25l . is a counterintuitive phenomenon occurring in a large variety of non-linear systems, 
whereby the addition of noise to a weak periodic signal causes it to become detectable or 
enhances the amount of transmitted information through the system |[m[T2l[T3l[T5l[T6l[r7l 
[ia[Ml2iinil23[33[Ml[3l[Ml[3a[3l][32l|2i|^ when SR occurs, the response of the 
system undergoes resonance-like behavior as a function of the noise level. In spite of the 
fact that initially this phenomenon was restricted to bistable systems, it is well known that 
SR appears in monostable, excitable, and non-dynamical systems. 

Here we report on experiments conducted on the response of A^. viridula individuals 
to sub-threshold signals. The non-monotonic behavior of SDM, as a function of the noise 
intensity (see Fig.|2j)), can be considered the hallmark of the threshold stochastic resonance 
(TSR). This phenomenon is well described by an extremely simple system, shown in Fig[3l 
and characterized by: (i) an energetic activation barrier (threshold); (ii) a weak coherent 
input such as a periodic signal (sub-threshold signal); (iii) a source of noise which is in- 
herent to the system, or is added externally to the deterministic input |[Tn[T2l[T3l . Since 
the three ingredients ai^e often present in nature and the idea of the existence of a thresh- 
old is quite intuitive, TSR has migrated into many different fields, so that during the last 
decades a considerable amount of literature on this subject has appeared in several areas 
of science and engineering. We have simulated a system with a threshold 0.045 V and a 
subthi^eshold signal of RMS amplitude 0.031 V (a. u.), obtained by the recorded female 
calling song. In Fig. |4] we show the output signal, and the corresponding PSD, for three 
different levels of noise added to the subthreshold deterministic signal (calling song). In the 
Figs. |4^, 11]:, 13^ we have rescaled the values of the signal amplitude in such a way that the 
zero value corresponds to the threshold value. For low noise intensities the signal crosses 
the threshold (dashed line in Fig. |4H) very rai^ely, and in the con^esponding PSD (Fig. |4j5) 
no frequency shows any significant power enhancement. By increasing the noise level the 
threshold crossings become more frequent (Fig. St) and the PSD appears to take a larger 
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threshold crossing spikes 



Figure 3. Evolution driven by a sinusoidal function plus noise, (a) Time series generated 
by consecutive pulses (dashed line: thi^eshold level, solid line: mean value of the periodic 
signal); (b) Temporal sequence of threshold crossing events. 
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Figure 4. Temporal evolution of the simulated calling signal over the threshold for noise 
intensity D = 2.6 • lO'^ V"^ (a), D = 1.30 • 10^^ (b) and L> = 1.0 • lO"^ V'^ (c). The 
corresponding power spectral densities are shown in panels b, d, f. The threshold level is 
s^i^ = 0.045 V and RMS amplitude of the subthreshold signal is 0.031 V"^. In the figures 
(a), (c) and (e) we have rescaled the values of the signal amplitude in such a way that the 
zero value corresponds to the threshold value. 



value for / = 102.5 Hz (Fig. |4jl), that is the dominant frequency contained in the input 
signal. A further increase of noise intensity produces a degradation of the signal, a loss of 
coherence in the temporal sequence (Fig. |4^) and a reduction of the main peak (/ = 102.5 
Hz) in the PSD (Fig.gf). The signal-to-noise ratio (SNR) at / = 102.5 Hz is reported in 
Fig. [5] For each value of the noise intensity we have performed 5000 numerical realizations. 
The noise intensity for which the SNR is maximum is Z) « 1.17 • 10"'^, which is very near 
the value of the noise intensity that maximizes the percentage of SDMs (see Fig. Wp)- The 
results obtained from our model suggest that in the biological system analyzed, stochastic 
resonance plays a key role, since it permits information to be extracted from a weak de- 
terministic signal, thanks to the constructive action of environmental noise. In other words 
there is a suitable noise intensity which maximizes the activating behavior of the green bugs 
and this effect can be described by the simplest threshold model which shows stochastic res- 
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4 LJ \ I \ I I I \ lJ \ I \ 

10"^ 10"'' 

Figure 5. Signal to noise ratio versus variance noise D of the output signal model when the 
input female calling song is subthi^eshold, at the dominant frequency / = 102.5. All the 
other parameters are the same of Fig. ID 



onance. In Fig.|5]we report also the best fitting curve of the simulations (cross points in the 
figure) obtained by the formula of the SNR for a single frequency coherent signal ifTTI 

SiVi? = clog[-^ exp(-A)], (1) 
where a = 6.6 • 10-^ b = 1.70 • 10"^ and c = 2.18. 



3. The Polymer Chain Model and MD Simulations 

The polymer is modeled by a semi-flexible linear chain of N beads connected by harmonic 
springs fTTl . Both excluded volume effect and van der Waals interactions between all beads 
are kept into account by introducing a Lennard- Jones (LJ) potential. In order to confer a 
suitable stiffness to the chain, a bending recoil torque is included in the model, with a rest 
angle = between two consecutive bonds. The total potential energy of the modeled 
chain molecule hU = L^Har + t^Bend + f^Lj with 

N-l 

f/Har = Yl Kr{r^,^+l " df (2) 
i=l 
N-l 

Usend = ^ Ko{9i_i^i+i — O^f (3) 
1=2 
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where is the elastic constant, rij the distance between particles i and j, d the equilibrium 
distance between adjacent monomers, Kg the bending modulus, elj the LJ energy depth and 
(T the monomer diameter. The effect of temperature fluctuations on the dynamics of a chain 
polymer escaping from a metastable state is studied in a two-dimensional environment. The 
polymer translocation is modeled as a stochastic process of diffusion in the presence of a 
potential barrier having the form: 

U^^t{x) = ax^ - bx^ (5) 

with parameters a = 3 • 10~^ and 6 = 2- 10~^, as already adopted in Ref. |[69l . A three- 
dimensional view of [/Ext is plotted in Fig. |6] The drift of the i"^ monomer of the chain 



Figure 6. 3D- view of the potential energy [/Ext> which is included in our system to simulate 
the presence of a barrier to be surmounted by the polymer. A sketch of the translocating 
chain molecule is shown. 



molecule is described by the following overdamped Langevin equations: 

— = — ^ ^ + Vl^Cx + A cos{ujt + (p) (6) 



dyi _ dUjj 
dt dy 



+ Vd^y (V) 



where Uij is the interaction potential between the i^^ and j^^ beads, and are white 
Gaussian noise modeling the temperature fluctuations, with the usual statistical properties, 
namely = and + r)) = 6(^^,1)^^) for (^,1 = x,y). A and uj are 

respectively the amplitude and the angular frequency of the forcing field and (p a randomly 
chosen initial phase. In our simulations, the time t is scaled with the friction parameter 7 as 
t = tr/j, where tr is the real time of the process. The standard Lennard- Jones time scale 
is tlj = (mcr^/eLj)^/^, where m is the mass of the monomer. A bead of a single-stranded 
DNA is formed approximately by three nucleotide bases and then cr ~ 1.5 nm and m 936 
amu |[65l . Orders of magnitude of the quantities involved in the process are nanometers for 
the characteristic lengths of the system (polymer and baiTier extension) and microseconds 
for the time domain. A set of 10'^ numerical simulations has been performed for different 
values of the frequency of the forcing field and two values of the noise intensity D, namely 
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D = 1.0,4.0. The values of the potential energy parameters are: = Kg = 10, eLj = 
0.1, fj = 3 and d = 5, in arbitrary units (AU). The amplitude of the electric forcing field 
is ^ = 5 • 10^^ in AU, because it is scaled with 7. The number of monomers is 20. 
The initial spatial distribution of the polymer is with all monomers at the same coordinate 
xq = 0, coiTcsponding to the local minimum of the potential energy of the barrier. Every 
simulation stops when the x coordinate of the center of mass of the chain reaches the final 
position at = 15. 
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Figure 7. MFTT vs. frequency of the forcing field for two different values of the noise 
intensity D. The values of the potential energy parameters ai^e: = Kq = 10, clj = 0.1, 
(7 = 3 and d = 5, in arbitrary units (AU). The amplitude of the electric forcing field is 
A = 5 • 10"^ (AU). The number of monomers N is 20. 




3.1. Results and discussion 

The MFTT shows three different translocation regimes as a function of the frequency 
(Fig. |7]). In the low frequency domain {uj < lO"'^), the period of the forcing field oscil- 
lations is very long with respect to the typical values of the mean crossing time of the chain 
molecule. In this regime the MFTT is equal to the average of the crossing times over up- 
per and lower configurations of the barrier, and the slowest process determines the value 
of the mean crossing time. As a consequence, the MFTT increases and we observe long 
tails in the probability density function (PDF) shown in Fig. [8^. In the high frequency 
domain {uj > 10^^), a saturation of the translocation time is obtained. In this case, very 
rapid oscillations act on the polymer motions as the mean potential, i. e. the static field, 
and therefore the MFTT becomes equal to that obtained without any additional periodic 
driving. In other words, the polymer chain feels the average potential barrier. For interme- 
diate frequencies (10~^ < lo < 10^^), the crossing event is strongly correlated with the 
potential oscillations and the MFTT vs. uj exhibits a minimum at a resonant oscillation rate. 
This frequency region coiTcsponds to periods of oscillations which are of the same order 
of magnitude of the mean time the polymer takes to cross a static barrier with its shape 
corresponding to the lowest configuration of the oscillating potential. In other words, the 
potential remains around its lowest configuration for enough time to allow the polymer to 
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Figure 8. Probability density function (PDF) of the first translocation time (FTT). Each 
panel shows two PDFs, each one chaiacterized by a specific value of the noise intensity. 
The three panels differs for the frequency of the forcing field. The panel (a) shows the time 
distribution in the low frequency region. The long tails indicate that the polymer crosses the 
potential barrier with a longer mean time. In the panel (b) the FTTs are distributed towards 
shorter values, because of the lowest time scale characterizing the translocation process in 
the resonant activation regime. The panel (c) shows the probability distribution for the high 
frequency domain, where the time scale is the same that is characterized by the presence of 
a static potential. 
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exit and, even in the case of an initially high or intermediate value of the height of the bar- 
rier, the potential feature turns into the lowest configuration within a sufficiently short time 
lag to facilitate the translocation process. The polymer, driven by a periodic field oscillating 
at a period comparable with a characteristic time of the crossing dynamics, reaches a res- 
onant regime that accelerates the translocation process. For each of the frequency values, 
the thermal noise intensity D is able to speed up or slow down the crossing process, as 
described by the three frequency regions (Fig. IT]) and the corresponding translocation dy- 
namics lISTIl . The probability density function of the first translocation time (FTT) is shown 
in Fig. [8] for three frequency values characterizing the different dynamical domains. Each 
panel shows two PDFs, each one characterized by a specific value of the noise intensity. 
In the resonant activation regime (Fig. [Sh) the PDFs do not present the long tail at higher 
crossing times, observed in Fig. [8^. Consequently, the MFTT reduces its value. The PDF 
assumes an interesting two-peaks structure that suggests the presence of two characteristic 
times of translocation. This feature, being present both at low and high noise intensity, can 
be ascribed to two different translocation dynamics of the polymer chain surmounting the 
barrier. In the high frequency domain (Fig. [8]:) the PDFs show the characteristic feature of 
the static potential case. 



4. Verhulst Model with Levy White Noise Excitation 

In considering how the population density x (t) may change with time t, Verhulst proposed 
the following equation 

dx ( x\ 

rx[\- - ]. (8) 



di \ Vl, 

where there is the Malthus term with the rate constant r and a saturation term with the 
factor, which is the upper limit for the population growth due to the availability of the 
resources. 

Really the parameters r and Vl are not constant. In fact the parameter r changes ran- 
domly due to season fluctuations, and the parameter fluctuates due to the environmental 
interaction which causes the random availability of resources. As a consequence we have 
the following stochastic Verhulst equation 



dx 

- = r(t)x 



1 



(9) 



In the context of macromolecular self-replication, the model equation with constant 
and a white Gaussian noise in r (t), was numerically studied in Ref. HI 02 II and the critical 
slowing down, i.e. a divergence of the relaxation time at some noise intensity, was found. 
Later Jackson and co-authors 111031 investigated the same model, by analog experiment 
and digital simulations. They analyzed specifically in detail the nonlineai^ relaxation time 
defined as ifTOil 



/,°°[(x(t))-(x(oo))]dt 



(10) 



x(0) - (x(oo)) 

and did not observe the critical slowing down. They explained this discrepancy by the in- 
correct approximate truncation of the asymptotic power series for T used in Ref. II102I . The 
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stability conditions were derived in Ref. II105I . Similar investigations for colored Gaussian 
noise r{t) were performed in Ref. 111061 . where a monotonic dependence of the relaxation 
time and the correlation time on the noise intensity was found. The generalization of Eq. Q, 
to study a Bernoulli-Malthus-Verhulst model driven by a multiplicative colored noise, was 
analyzed recently in Ref. II107I . 

The evolution of the mean value in the case of Eq. Q with constant r and white Gaus- 
sian noise excitation /3 (t) = r/n (t) was considered in Refs. ll85l [T08l[T09l[TT0l[TTTl[TT2l . 
In Refs. IIIIIIIII2I the authors, using perturbation technique, obtained the exact expansion 
in power series on noise intensity of all the moments and found the long-time decay of 
In Ref. f85| the authors derived the long-time behavior of all the moments of the 
population density by means of an exact asymptotic expansion of the time averaged process 
generating function, and found the same asymptotic behavior of at the critical point. 
This very slow relaxation of the moments near the critical point is the phenomenon of crit- 
ical slowing down. 

In the present chapter, using the previously obtained results for a generalized Langevin 
equation with a Levy noise source II113II1141 . we investigate the transient dynamics of the 
stochastic Verhulst model with a fluctuating growth rate and a constant value for the sat- 
uration population density Q, that is = 1. The exact results for the mean value of the 
population density and its nonstationary probability distribution for different types of white 
non-Gaussian excitation r (t) are obtained. We find the interesting noise-induced transi- 
tions for the probability distribution of the population density and the relaxation dynamics 
of its mean value for Cauchy stable noise. Finally we obtain a nonmonotonic behavior of 
the nonlinear relaxation time as a function of the Cauchy noise intensity. 



5. Stochastic Verhulst Equation with Non-Gaussian Fluctua- 
tions of Growth Rate 

Let us consider Eq. ^ with a constant saturation value Q = 1, namely 

dx 

— =r{t)x{l-x). (11) 
After changing variable y = ln[x/(l — x)], we obtain 

y{t) = y{0)+ [ r (r) dr 
Jo 

and the exact solution of Eq. (fTTI) is 

x{t)=(l + ^—^eKpl- [\{T)dT\) , (12) 



I Jo 

where xq = x (0) . Now by substituting in Eq. ([121) the following expression for the random 
rate r{t) 

r (t) = r + e (t) , (13) 
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where r > and ^ (t) is an arbitrary white non-Gaussian noise with zero mean, we can 
rewrite the solution ([T2l ) as 

^(t)= (^l + lz^e-*-M*)^"\ (14) 

Here L (t) denotes the so-called Levy random process with L (0) = 0, and ^ (t) = L {t). 
As it was shown in Refs. Ill 13111 14111151 . Levy processes having stationaiy and statistically 
independent increments on non-overlapping time intervals belongs to the class of stochas- 
tic processes with infinitely divisible distributions. As a consequence, the characteristic 
function of L (t) can be represented in the following form (see Eq. (6) in 111131 ) 

(e»«'>) = exp {i £ , (.) a.} . (15) 

where p {z) is some non-negative kernel function. The case p (z) = 2D6 (z) corresponds 
to a white Gaussian noise excitation ^ (t), while for a symmetric Levy stable noise (t) 
with index a we have a power-law kernel p{z) = Q |^;|^~°, with < a < 2. 

In the model under consideration the stationary probability distribution has; (i) a singu- 
larity at the stable point x = 1 for white Gaussian noise; and (ii) two singularities at both 
stable points a; = and x = 1 for Levy noise. To analyze the time behavior of the probabil- 
ity distribution in the transient dynamics it is better not to use the Kolmogorov equation for 
the probability density P {x, t), but rather the exact solution (fT4l) . Using the standard the- 
orem of the probability theory regarding a nonlineai" transformation of a random variable, 
we find from Eq. ([T4l) 

1 / rn _ Tr.^ \ 

rt,t], (16) 



Pix,t) ^ 



X (1 — x) 





"(1 - xo)x' 


Xq (1 - X)_ 





where Pl {z,t) is the probability density con^esponding to the characteristic function ([T5] ). 
For a white Gaussian noise ^ (t), this distribution reads 



PL{z,t) = -^exp{-^\. (17) 



I 4L>tj' 

The time evolution of the probabiUty distribution P {x, t) for D = 0.3, r = 2, and xq = 0.1 
is plotted in Fig. |9l 

As it is easily seen, the maximum of the unimodal distribution with initial position at 
X = 0.1 shifts with time towards the stable point at x = 1. At the same time, as it follows 
from Eqs. (fT6l) and (fTTl ). for all t > we have 

lim P{x,t)= lim P{x,t) = 0. (18) 

The same picture is observed for another kernel function p{z) = Kz/ (2sinhz) 
{K > 0) , corresponding to a Levy process ry {t) with finite moments and the following 
probability density of increments 



2Kt-i fKt iz\^fKt 



iz 
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Figure 9. Time evolution of the probability distribution of the population density for white 
Gaussian noise excitation with intensity D. The values of the parameters are: xq = 0.1, 
r = 2,D = 0.3. 



P{x,t) 




0.2 0.4 0.6 0.8 1 
Figure 10. Time evolution of the probability distribution of the population density in the 
case of Levy noise with distribution ( [T9l ). The values of the parameters are: xq = 0.1, 
r = 2,K = 0.2. 
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where T (x) is the Gamma function. The corresponding time evolution of the probability 
distribution P {x, t) for K = 0.2, r = 2, and xq = 0.1 is shown in Fig.fTOl 

A different situation we have for a Cauchy stable noise ^ (t) with constant kernel 
p{z) = Q {a = 1). After evaluation of the integral in Eq. ([T5] ). the probability density 
of the Levy process increments takes the form of the well-known Cauchy distribution II115II 

PL{z,t)= (20) 



where Di = ttQ is the noise intensity parameter. In such a case from Eqs. (fT6b and (1201 ) for 
alH > we find 

oo. (21) 



lim P (x, t) 

r-s>0+ 



lim P (x, t) 



As a result, from an initial delta function we immediately obtain a trimodal distribution 
for i > and then after some transition time tc a bimodal one with two singularities at 
the stable points x = and x = 1 (see Figs. [TTI - [T3] ). We should note that the transition 



P{x,t) 




Figure 1 1 . Time evolution of the probability distribution of the population density in the 
case of white Cauchy noise excitation. The values of the parameters are: xq = O.l, r = 2, 
Di = 0.7. 

from trimodal to bimodal distribution is a general feature of the model in the presence of a 
Cauchy stable noise, and it is not limited to some range of parameters. In fact, from Eq. (|2T]) 
and a delta function initial distribution inside the interval (0,1), this transition always takes 
place. 

In the following Figs. [12] and [13] we show the time evolution of the probability distribu- 
tion of the population density for two other values of the noise intensity, namely Di = 1.2 
and Di = 1.7. As the noise intensity increases the probability distribution shows two 
singularities near x = and x = 1 with different amplitude. 

This transition in the shape of the probability distribution of the population density is 
due to both the multiplicative noise and the Levy noise source. Using Eqs. ([T6] ) and (|20] ) 
and equating to zero the derivative of P{x, t) with respect to x, we obtain the following 
condition for the extrema in the range < x < 1, and particularly for a minimum in the 
same interval 

z(x.t) 1 
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Figure 12. Time evolution of the probability distribution of the population density in the 
case of white Cauchy noise. The values of the parameters are xq = 0.1, r = 2, Di = 1.2. 




Figure 13. Time evolution of the probability distribution of the population density in the 
case of white Cauchy noise. The values of the parameters are xq = 0.1, r = 2, Di = 1.7. 
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with 



z{x, t) = In 



' (1 - xo)x ' 
xq (1 - x)_ 



rt . 



(23) 



This condition can be solved graphically by finding the intersection between the functions 
Ui = z{x, t) I (z(x, t)^ + (D\tY) and y2 = x — 1/2. This is done in the following Figs. [14]- 
[161 where the function yi is plotted for thi^ee different values of time and noise intensity. In 
each figure the black blue curve (color on line) corresponds to the critical value of time 
for which we have a noise induced transition of the probability distribution of the popula- 
tion density from trimodal to bimodal, that is from two minima and one maximum to one 
minimum inside the interval < a; < 1. The appearance of one minimum in the probability 
distribution is the signature of this tiansition. 




Figure 14. Plots of both sides of Eq. ([22l) (white Cauchy noise): function yi (sohd curves), 
function 7/2 (dashed curve), for three values of time, namely: t = 0.8, 1.75, 2.5. The critical 
time is tc = 1.75 (black blue curve). The values of the other parameters are: xq = 0.1, 
r = 2,Di = 0.7. 




Figure 15. Plots of both sides of Eq. ([22] ) (white Cauchy noise): function yi (solid curves), 
function y2 (dashed curve), for three values of time, namely: t = 0.9, 1.3, 3. The critical 
time is tc = 1.3 (black blue curve). The values of the other parameters are: xq = 0.1, 
r = 2,Di = 1.2. 

The three values of the critical time tc corresponding to the three values of the Levy 
noise intensity investigated are: Di = 0.7, tc = 1.75; D2 = 1.2, tc = 1.3; = 
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1.7, tc = 0.95. One rough evaluation of the critical time tc is obtained by putting equal 
to 1 the scale parameter of the Cauchy distribution of Eq. (l20l ). that is Tc ~ l/Di. The 
critical time tc is the time at which the maximum and one minimum of the probability 
distribution (see Figs. [TTI - [T3] ) coalesce in one inflection point and in this point x the function 
1/2 = X — 1/2 becomes tangent at the function yi (see Figs. [T4] - [T6l ). It is interesting to note 
that the critical time tc decreases with the noise intensity Di . This is because by increasing 
the noise intensity, more quickly the population density reaches the two points near the 
boundaries x = and x = 1. 



1 




-1 



Figure 16. Plots of both sides of Eq. (l22l ) (white Cauchy noise): function yi (solid curves), 
function y2 (dashed curve), for three values of time, namely: t = 0.4, 0.95, 1.5. The critical 
time is tc = 0.95 (black blue curve). The values of the other parameters are: xq = 0.1, 
r = 2,Di = 1.7. 



6. Nonlinear Relaxation Time of the Mean Population Density 

It must be emphasized that to find the time evolution of the mean population density one 
can use two different approaches. The first one was proposed in Ref. 111031 . According to 
the exact solution ([141 ) of the Verhulst equation (ITTI ). we can rewrite this expression in the 
following form 

x(t) = /(e-'-*"^W), (24) 

where 

.f(q)=(l + - 

Xq 

Then, by expanding the smooth function (1251 ) in a standai'd Taylor power series in q around 
the point g = we have 

^-^ n! 

After substitution of Eq. ( |26l ) in (l24l ) and averaging we obtain 

n=0 



f{q)=(l + ^^q{t)] . (25) 
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or, in accordance with Eq. (fTSl ). 



{xit)) = E „, (28) 



n=0 



r /■+°° e""^ - 1 + nsinz , , . 
X exp|t J 2 p[z)dz^. 



For white Gaussian noise ^ (t) with kernel p (z) = 2D5{z) we obtain from Eq. (1281 ) the 
following asymptotic series 



(^(t)) = yZ!lWe^*-^— (29) 



n=0 



By considering a finite number of terms in this expansion leads to a wrong conclusion about 
the critical slowing down phenomenon in such a system, as found in Ref. II 1021 . The exact 
result is obtained, of course, by summing all the terms in Eq. (|29l ). Moreover, for most of 
the kernels p{z) the integral in Eq. ( |28] ) diverges. Thus, this approach is inappropriate for 
our purposes, and it is better to use the direct averaging in Eq. ([T4l) . Therefore, using this 
second approach we have 

(x (t)) = (l + e-'-*--') ' Pl {z, t) dz. (30) 

Let us consider now different models of white non-Gaussian noise ^ {t). We start with 
the white shot noise 

i{t) = Y,'^i5{t-U) (31) 

i 

having the symmetric dichotomous distribution of the pulse amplitude P{a) = 
[5 (a — ao) + 5{a + oq)] /2, mean frequency v of pulse ti^ain, and kernel p {z) = i/z^P {z). 
From Eq. ( [T5] ) we have 

' ^iuL(t)\ _ ^-ut(l-cosaou) ^^2) 



By making the reverse Fourier transform in Eq. (1321 ) we find the probability distribution of 
the corresponding Levy process 

+00 

PL{z,t) = e-'' Yl In(.iyt)6{z-nao), (33) 

71= — 00 

where /„ (x) is the n-order modified Bessel function of the first kind. The relaxation of the 
mean population density {x{t)) is shown in Fig. [T7] According to the Eqs. (l30l ) and (1331 ) 
the stationary value of the population density in such a case is {x)^^ = 1, but the relaxation 
time ([Tol l increases with increasing the mean frequency of pulses. 

For white non-Gaussian noise with the kernel p{z) = Kz/ (2 sinh z) we observe a sim- 
ilai" transient dynamics, which is shown in Fig.[T8l We have the same stationary value (x)^^, 
and the relaxation time T increases with increasing the parameter K, which is proportional 
to the noise intensity. 
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Figure 17. Nonlinear relaxation of the mean population density in the case of white shot 
noise excitation, for three values of the mean frequency u, namely u = 0.3, 3. The values 
of the other parameters are: xq = 0.1, r = 2, ao = 1. 



Finally, in the case of white Cauchy noise ^ {t) we obtain interesting exact analytical 
results. First of all, substituting Eq. ( [201 ) in ( [30l ) and changing the variable z = Dity under 
the integral, we obtain 

{x (t)) = i \l + l^e-*(^+^^^)l . (34) 



For the stationary mean value {x)^^ we find from Eq. (l34l ) 

where 1 {x) is the step function. After evaluation of the integral in Eq. (1351 ) we obtain finally 

,,11 r 

\^)st = o + - arctan — . (36) 
z vr JJi 



As it is seen from Fig. [191 and Eq. (I361 ). for small noise intensity Di, with respect to the 
value of the rate parameter r = 2, the stationary mean value of the population density is 
approximately 1, as for the other white non-Gaussian noise excitations considered. But for 
large values of Di, this asymptotic value, which is independent from the initial value of 
population density xq, tends to 0.5. 

It is interesting also to analyze, for this case of white Cauchy noise, the dependence 
of the relaxation time T from the noise intensity Di. Substituting Eq. ([34l ) in (ITOt and 
changing the order of integration, for initial condition xq = 0.5, we ai^e able to calculate 
analytically the double integral in t and in y obtaining the final result 

vrln2 

(37) 



■l + Z)f/r2)arccot(L>i/r) 



We find a nonmonotonic behavior of the relaxation time T versus the noise intensity Di 
with a maximum at the noise intensity Di = 0.75, as shown in Fig.l20l This nonmonotonic 
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Figure 18. Nonlinear relaxation of the mean population density in the case of Levy noise 
with distribution (l20l ). for three values of the parameter K, namely K = 0.4, 2. The values 
of the other parameters are: xq = 0.1, r = 2. 

behavior is also visible for another initial position xq = 0.1 in Fig. [19] Here the relaxation 
time to reach the stationary value of mean population density (x) increases from very low 
noise intensity {Di = 0.2) to moderate low intensity {Di = 2), while decreases for higher 
noise intensities {Di = 7). This is also due to the dependence of (x)^^ from the noise 
intensity Di (see Eq. (|36l)). We note that this nonmonotonic behavior of the relaxation time 
T is related to the peculiarities of the transient dynamics of the mean population density 
and it will be object of further investigations. 

7. Conclusions 

In this contribution we have studied the effects of random fluctuations, i.e. noise, in the 
vibrational communications occurring during the mating of A^. viridula, i.e. the green bug. 
In our experimental work we analyzed the behavioural response of different individuals of 
A^. viridula to a deterministic signal (calling song), measuring for these individuals the the 
threshold of the neural activation. Afterwards, we analyzed the green bug response when 
a sub-threshold deterministic signal is added with an external noise source. By using the 
Source-Direction Movement ratio as indicator of positive response to the external signal, we 
observe that the behavioural activation of the insects is characterized by a non-monotonic 
behaviour as a function of the noise intensity D, with a maximum at D = D^pt « 1.30 • 

maximizes the efficiency of the sexual communication among 
individuals of Nezara viridula (L), and therefore represents the optimal noise intensity 
during the mating behaviour of these insects. The non-monotonic behaviour observed in 
the insect response as a function of the noise intensity is the signature of the thi^eshold 
stochastic resonance (TSR) lITSl . By using a threshold model we obtained numerical results 
for the threshold crossing, which coiTcsponds in our model to the behavioural activation, 
finding a theoretical value for the optimal noise intensity. Experimental and numerical 
results are compaied, finding a good agreement between the values of the optimal noise 
intensity obtained by the experimental work and model (see Figs.|2tb) and[5]l. 



26 



B. Spagnolo, D. Valenti, S. Spezia et al. 



J <x(0> 




2 4 6 8 10 



Figure 19. Nonlinear relaxation of the mean population density in the case of white Cauchy 
noise, for three values of the noise intensity Di, namely Di = 0.2, 2, 7. The values of the 
other parameters are: = 0.1, r = 2. 
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Figure 20. Nonmonotonic behavior of the nonlinear relaxation time T as a function of the 
white Cauchy noise intensity Di. The values of the other parameters are: xq = 0.5, r = 2. 

We analyzed the influence of an external oscillating driving field on the translocation 
dynamics of short polymers embedded in a noisy environment. We simulate the transloca- 
tion process by letting the polymer to cross a potential barrier starting from a metastable 
state, in the presence of thermal fluctuations. The mean translocation time as a function of 
the frequency of the driving force shows a nonmonotonic behaviour, with the noise intensity 
acting as a scaling factor of the values of the crossing times. The forcing periodic electric 
field jointly with the temperature of the system can be able to speed up or slow down the 
polymer translocation. In this view, the oscillating electric field constitutes a tuning mecha- 
nism to select a suitable translocation time of the polymer. This feature may have important 
biological effects on the cell metabolism, for example, during a cancer targeted therapy. 

Finally, the transient dynamics of the Verhulst model, perturbed by arbitrary non- 
Gaussian white noise, has been investigated. This well-known equation is an appropriate 
ecological and biological model to describe closed-population dynamics, self -replication of 
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macromolecules under constraint, cancer growth, spread of viral epidemics, etc... By using 
the properties of the infinitely divisible distribution of the generalized Wiener process, we 
analyzed the effect of different non-Gaussian white sources on the nonlinear relaxation of 
the mean population density and on the time evolution of the probability distribution of the 
population density. We obtain exact results for the nonstationary probability distribution in 
all cases investigated and for the Cauchy stable noise we derive the exact analytical expres- 
sion of the nonlinear relaxation time. Due to the presence of a Levy multiplicative noise, 
the probability distribution of the population density exhibits a transition from a trimodal 
to a bimodal distribution in asymptotics. This transition, characterized by the appearance 
of a minimum, happens at a critical time tc, which can be roughly evaluated as tc ~ l/^i 
(where Di is the noise intensity) and exactly evaluated from the condition (l22l) . Finally 
a nonmonotonic behavior of the nonlinear relaxation time of the population density as a 
function of the Cauchy noise intensity was found. 
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